Hybrid reciprocal space for X-ray diffraction in epitaxic layers 
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Even after several decades of systematic usage of X-ray diffraction as one of the major analyt- 
ical tool for epitaxic layers, the vision of the reciprocal space of these materials is still a simple 
superposition of two reciprocal lattices, one from the substrate and another from the layer. In this 
work, the general theory accounting for hybrid reflections in the reciprocal space of layer/substrate 
systems is presented. It allows insight into the non-trivial geometry of such reciprocal space as 
well as into many of its interesting properties. Such properties can be further exploited even on 
conventional-source X-ray diffractometers, leading to alternative, very detailed, and comprehensive 
analysis of such materials. 



I. INTRODUCTION 

The capability of growing thin layers of single-crystals 
onto one face of another single crystal has made possible 
many fundamental achievements in semiconductor tech- 
nology. Epitaxic growth is today one of the most impor- 
tant and basic processes used in manufacturing nanos- 
tructured devices. Multilayered materials, such as su- 
perlattices and quantum wells, or even quantum wires 
and dots, require epitaxy at some stage of their prepara- 
tion procedures. X-ray diffraction has been the primary 
tool for structural analysis of epitaxic layers, with the as- 
sociated techniques and machinery following closely the 
needs of the semiconductor industry. 

A quarter of a century ago, when using a diver- 
gent X-ray source and photographic films to record the 
layer/substrate diffraction lines [the simplest possible 
setup (Chang, 1980) to measure lattice mismatch of epi- 
layers], Isherwood et al. (1981) reported the observation 
of extra features, a kind of short line, appearing all over 
the recorded images. Such features were sequences of 
consecutive Bragg reflections in both single-crystal lat- 
tices, and were named hybrid reflections. Latter, the 
phenomenon was quantitatively described and methods 
to exploit its properties were suggested (Morelhao & Car- 
doso, 1991; 1993a; 1993b). However, to probe the excite- 
ment conditions of hybrid reflections precisely, coUimated 
X-ray beam setups would be necessary, such as those 
commonly found in most synchrotron facilities where the 
beam can be highly coUimated in two orthogonal direc- 
tions (Morelhao et al, 1991, 1998, 2002, 2003). Such 
requirements have created a huge barrier for the system- 
atic usage of this peculiar phenomenon in the technology 
of semiconductor devices. 

Even after several decades of using X-ray diffraction as 
one of the major analytical tools for epitaxic layers, the 
vision of the reciprocal space of these materials is still a 
simple superposition of two reciprocal lattices: one from 
the epilayer and the other from the substrate. Diffraction 
conditions generating any other extra feature have been 
avoided since they could not be explained within this sim- 
plistic vision of the reciprocal space. The incident-beam 
optics available for conventional X-ray sources have, in 



the past, seemed inappropriate (considering the require- 
ment for a highly coUimated beam in two orthogonal 
directions) to investigate azimuthal-dependet features; 
consequently, the analysis of epilayers by standard X-ray 
diffraction techniques would be compromised when hy- 
brid reflections are excited. This work is an attempt to 
change this scenario. Introducing a reciprocal-space de- 
scription of hybrid reflections opens the possibility of ex- 
ploiting in detail the properties of this phenomenon with- 
out the necessity of using synchrotron facilities. In other 
words, instead of avoiding hybrid features when using 
tube-source diffractometers, exciting them via standard 
reciprocal-space probing techniques can lead to alterna- 
tive methods for analyzing epilayered materials. Here, 
besides presenting the theory and discussing some prop- 
erties of the hybrid reciprocal lattice, experimental ex- 
amples are given regarding the type of information that 
can be accessed by analyzing the phenomenon properly 
on commercial diffractometers. 



II. HYBRID RECIPROCAL LATTICE THEORY 

Any three-dimensional reciprocal lattice gives rise to a 
phenomenon known as rt-beam diffraction (Colella, 1974; 
Chang, 1984; Weckert & Hiimmer, 1997). AUhough it 
can change the relative strength of Bragg reffections, 
no extra features are generated in the reciprocal space 
since sums of diffraction vectors always end up at a 
reciprocal-lattice point (RLP). On the other hand, when 
two distinct reciprocal lattices are superposed, as in epi- 
layer/substrate systems, sum of diffraction vectors may 
end up at an empty position of the reciprocal space. This 
occurs when one diffraction vector in the sum does not 
belong to the same lattice as the others. In this case, hy- 
brid reciprocal-lattice points (HRLPs) are generated, as 
illustrated in Fig. 1 and described mathematically below. 

Three-beam X-ray diffractions in single crystals are ex- 
cited when the incident beam (wavevector k and wave- 
length A), fulflUs two Bragg conditions: 
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FIG. 1: In epilayer/substrate systems, superposition of the 
substrate (circles) and epilayer (ellipses) reciprocal lattices 
gives rise to inter-lattice rescattering processes, as illustrated 
by the hybrid reciprocal-lattice point (HRLP) with diffrac- 
tion vector P' = Ms + Nl- To excite one chosen HRLP, the 
crystal must be at specific azimuthal ip angles as detailed in 
Appendix A. 

Since P = M + N, we also have 

k- N = -N ■ N/2- N ■ M (3) 

where P, M, and TV are diffraction vectors of the pri- 
mary, secondary and coupUng reflections, respectively. 
The primary reflection is the one whose intensity is been 
monitored while the secondary reflection is brought to 
diffraction condition by the crystal azimuthal rotation 
around P, as in X-ray Renninger scanning (Renninger, 
1937). Other cases of n-beam diffractions with n > 3 will 
be treated here as coincidental three-beam diffractions. 

Only equations H]) and ^ are in fact necessary to 
predict three-beam diffractions in a single lattice (Cole 
et ai, 1962; Caticha, 1969), which can be either the epi- 
layer or the substrate one. However, there are several 
other similar diffraction processes: the above-mentioned 
hybrid reffections, whose secondary and coupling reflec- 
tions do not belong to the same lattice. To predict what 
should be the exact incident-beam direction for excit- 
ing one of such inter-lattice rescattering process, equa- 
tions ([2]) and ^ are more suitable to this purpose, as 
demonstrated elsewhere for the case of satellite reflec- 
tions (Morelhao et al, 2003) and summarized here in 
Appendix A for the case of epilayer/substrate systems. 
Accounting for all possible rescatterings leads to a re- 
ciprocal space that is highly populated with HRLPs and 
much more complex than that obtained by just superpos- 
ing both epilayer and substrate reciprocal lattices. This 
hybrid reciprocal space has been neglected; its features 
remain unexplored; all knowledge on this matter is found 
in reports of a few accidental observations (Hayashi et 



al. 1997, Domagala et al, 2006) sometimes investigated 
(Morelhao et ai, 2003), but most of the time avoided in 
order not to compromise the system characterization by 
standard diffraction techniques, such as rocking curves 
and reciprocal-space mapping in triple-axis goniometry. 

To visualize the hybrid reciprocal space, let us label 
the diffraction vectors of both lattices as 

Ms,L ha*gj^ + kb*gj^ + Icgj^ and 

Nl,s = {H - h)als + - k)bls + (L - £)cls 

where S and L subscripts stand for the substrate and 
epilayer reciprocal-lattice vectors, respectively, h, k, and 
£ are the Miller indexes of the secondary reflection, and 
the complete hybrid reciprocal space around one chosen 
HKL primary reflection of the substrate lattice, whose 
diffraction vector is P = Ha*g + Kbg + Lcg, only ap- 
pears by rotating the sample around P by 360°. The 
position of all features regarding the P vector is then 
given by 

AP = P' P = h'Aa* + k'Ab* + i'Ac* , (4) 

which is a sub-reciprocal-lattice of points with periodicity 
Ag* ~ g1 ^ 9s' for 9* — ^* ^ ^'^d c*, since h\ k' , and 
£' are integer numbers. They stand for either coupling- 
reflection {H — h, K — k, L — i) or secondary-reflection 
{h, k, £) indexes according to P' = Ms + Ni or P' = 
Ml + Ns, respectively. 

Although equation (|4]) gives the general aspect of the 
hybrid reciprocal lattice, there are a few restrictions 
that should be considered for each particular system. 
One is the direction of the secondary beam (wavevector 
kM = Ms.L + fe) that must cross the epilayer/substrate 
interface in order for its respective hybrid diffraction vec- 
tor, P' — Ms^L + ^L,S: to be measurable. In other 
words, if n is the interface normal direction pointing up- 
wards into the epilayer and ^ = n • kM/\kM\, we have 
that 

(h' k' n - / iH-h,K-k,L-e) ff e>o 

\ (Kkj) ff e<o. ^ ' 

Hence h'k'i' are the indexes of the Bragg reflection that 
occurs in the epilayer lattice. Refraction and total re- 
flection of the secondary beam, fc^f, at the interface is 
also a possibility to be taken into account mainly when 
^ w 0. In general, epilayers of quaternary alloys grown 
on miscut substrates present a relative tilt between both 
lattices, and this tilt must be considered when calculating 
Ag*. 

A. Properties of the hybrid reciprocal lattice 

One of the most interesting properties of the hybrid 
reciprocal lattice is that the relative positions AP of the 
HRLPs depend exclusively on the lattice mismatch be- 
tween both real lattices, i.e. their positions do not depend 
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on the X-ray wavelength. What change with A is the az- 
imuthal angle, (/?, where each HRLP is excited, as can be 
calculated by solving equations ^ and ^ ; see Appendix 
A. This implies that HRLPs are more easily excited with 
an X-ray beam that is poorly collimated in the axial di- 
rection, i.e. the direction perpendicular to the primary 
incidence plane: the plane that contains the substrate 
vector P, the X-ray source and the detector system. In 
commercial high-resolution diffractometers, the beam is 
highly conditioned in the incidence plane to about a few 
arc seconds, while the axial divergence is of the order of 
a few degrees (« 2°). On the other hand, in synchrotron 
facilities for X-ray diffraction, where the beam is well con- 
ditioned in both directions, the azimuthal positioning of 
the sample has to be more accurate if hybrid features are 
to be measured (Morelhao et ai, 2002; 2003). 

In Fig. 2, the general properties of hybrid reciprocal 
lattices are depicted. Just a few points aligned along 
the growth direction in the case of fully strained lay- 
ers (Fig. 2a), or a well defined three-dimensional lattice 
of points around the substrate reciprocal-lattice points 
occurs in the case of relaxed layers (Fig. 2b). Conse- 
quently, strain gradients along the layer thickness give 
rise to a hybrid lattice of rods instead of points, as 
illustrated in Figs. 2(c) and 2(d). Besides the fact 
that HRLPs are excited only at certain azimuthal po- 
sitions, their experimental observation via conventional 
reciprocal-space mapping techniques also require detec- 
tion optics with some angular acceptance in the axial 
direction. When a given HRLP is excited, its diffraction 
vector P' — AP -I- i-" is not on the primary incidence 
plane and neither is its diffracted beam k' — P' + k. 
Therefore, visible HRLPs are those whose diffracted 
beam, k', falls within the angular range of axial accep- 
tance of the detection system. Moreover, HRLP posi- 
tions on reciprocal-space maps correspond to projections 
of Ai-" on the incidence plane. In terms of longitudinal 
and transversal components of the maps, Qz and Qxy 
respectively, the HRLPs are seen at 



a) 



Qz — AP • z and Q^y = AP • k^ 



(6) 



where z = P/|P| and k^y — [k 
It is then possible to calculate 



{k- z)z]/\k~ {k- z)z\. 



S = \Qxytstn{Aip) 



(7) 



as the takeoff angle of k' from the incidence plane of the 
primary reflection. Aip = ip' — ip, and ip' is the azimuthal 
position of P' around P with respect to the same refer- 
ence direction for ip and in the same sense of rotation. 



1. Strain grading 

A unique property of hybrid reflections arises from the 
layer reflections since they are always Laue reflections, 
i.e. transmitted-diffracted beams where diffraction oc- 
curs through the entire layer thickness. In any other 
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FIG. 2: Properties of hybrid reciprocal lattice in epi- 
layer/substrate (001) cubic systems, (a) HRLPs near sym- 
metrical, 002, 004, and 006, and asymmetrical, 224, substrate 
reflections. In fully strained layers, Aa* = Ab* = 0, equa- 
tion (|4|, and hence the HRLPs are aligned along the (001) 
growth direction and they are distinguished only by their £' 
index, shown beside of each one of them (dark spots). The 
HRLP with £' = coincides with the substrate RLP marked 
by S. L stands for the layer RLP. (b) Relaxed layers where 
Aa* = A6* 7^ provide a three-dimensional hybrid lattice 
around either symmetrical, 004, or asymmetrical, 224, pri- 
mary substrate reflections, (c), (d) Strain grading along the 
layer thickness gives rise to elongated HRLPs towards the sub- 
strate one. HRLPs with index i' = lay on the layer in-plane 
direction, and (d) they are distinguishable from the substrate 
RLP if some relaxation occurs at the layer/substrate inter- 
face, the h'k' indexes are given. In the insets, a|| and a± 
stand for the in-plane (parallel) and out-plane (perpendicu- 
lar) unit cell parameters of the layer, respectively, ao is the 
substrate lattice parameter. In these examples, a± > ao, and 
only HRLPs in which j^| > 0.008, equation (O, are shown. 



X-ray diffraction technique, structural analyses of lay- 
ers are carried out by means of Bragg reflections: those 
with reflected-diffracted beams. Therefore, in principle, 
HRLPs are equally sensitive to lattice mismatch at both 
the top and the bottom of the epilayer. Hence, strain 
variation across the layer thickness should generate elon- 
gated HRLPs, i.e. HRL rods, oriented at specific direc- 
tions on reciprocal-space maps. Except for the particular 
cases discussed in Figs. 2(c) and 2(d), which are sensitive 
only to the in-plane strain grading, the general orienta- 
tions of the rods are affected by the elastic properties of 
the layer compound; for example, in the case of (001) 
growth on cubic systems with parallel, Ae||, and per- 
pendicular, Ae^, strain gradings. The h'k £' HRL rod 
around a symmetric primary reflection will be oriented 
at an angle XG where 
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Ae, 



AQz ^ -(.' 

AQxy ^ + fc'2 cos(A(^) A£|| 



(8) 
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FIG. 3: Reciprocal-space maps, 002 reflection, ZnSe (1/im thick) relaxed layer on GaAs (001), recorded with C\iKa\ radiation 
at different azimuthal ip angles, (a) (/p = 0, only the 002 substrate (S) and layer (L) RLPs are seen; the former is at Qz ~ 
2/ao = 0.35377A~^. (b) ip = 27.6°, (c) ip = 30°, and (d) ip — 45°. [110] is the in-plane reference direction for the (p rotation 
of the sample around the [001] direction. Visible HRLPs, in (b), (c), and (d), are numbered and identified by their transversal 
components Q^y in Table 1. White bars represent 1.6 x 10"^A"\ Mesh size: 2.48 x 10~^A-1 in Q^y per 3.27 x 10"^A-1 in 
Qz. 



assuming a tetragonal distortion of the layer unit cell 
varying as a function of the strain according to a = & = 
0^(1 +e||) and c = 0^(1 + e±); ol is the unstressed layer 
lattice parameter and, for isotropic materials, 



Aeii 



2iy 



1 



(9) 



where ly is the Poisson ratio. Being able to predict HRL 
rod orientation due to grading is important since there 
are other causes for the elongated shapes of the HRLPs, 
as better explained below. 



2. Anisotropic mosaicity 

Besides strain grading and finite layer thickness effects, 
the HRLPs may also present elliptical shapes due to mo- 
saicity in the epilayer. Most diffraction vectors taking 
part in hybrid reflections have both in-plane and out- 
plane components, and therefore, they are sensitive to 
the spatial misorientation of the crystallites, or mosaic 
blocks, building up the epilayer. In-plane rotation of the 
crystallites around the growth direction gives raise to a 
mosaicity of width rj\\ , while crystallite rotations around 
in-plane axes give rise to an out-plane mosaic width 'qji_. 
It is possible to compute effects of mosaicity on the shape 
and orientation of the HRLPs by projecting the trajec- 
tory of the reciprocal vector 

Aif — 77|| {z x r) sin tp + ri± [{z x r) x r] cos ip (10) 

on the incidence plane of the reciprocal-space maps as 
ip varies through 360°. f = H/\H\ and H is the layer 
diffraction vector in P' = Ms^l + -^l.s, i-e. either H = 



Ml ot H = Nl. The orientation angle xm is with 
respect to the direction of largest projection so that 



tan Xm 



AH z 
AH 



(11) 



when the value of AQ'j, + AQ^ is a maximum. 



III. RESULTS AND DISCUSSIONS 

Occurrence of experimental HRLPs is demonstrated 
in Fig. 3, obtained in a sample with a single l/um thick 
ZnSe epilayer on GaAs (001) substrate. Proper identifi- 
cation of the HRLPs is given in Table 1. The reciprocal- 
space maps were collected on a Philips X'Pert-MRD high 
resolution diffractometer: Cu tube. X-ray mirror, four- 
crystal asymmetric 220 Ge monochromator and three- 
bounce 220 Ge analyzer crystal; nominal spectral width 
AA/A = 5 x 10-5; X-ray beam divergences 0.005° and 2° 
in the incidence plane and in the axial direction, respec- 
tively. 

Direct and complete strain analysis of the epilayer is 
possible even from reciprocal-space maps of the sym- 
metric 002 GaAs reflection when HRLPs are excited. 
The Qxy components of the HRLPs in Figs. 3(b), 3(c), 
and 3(d) are well reproduced by using Aa* = Ah* — 
—5.15(12) X lO-'^A"-'^, and their longitudinal separation 
(an integer fraction of the layer-substrate reciprocal- 
lattice point distance) provides Ac* — —4.43(33) x 
lO-^A-i. Hence, assuming a Poisson ratio of = 0.5, 
which implies e_L/£|| — —2, the fully relaxed lattice pa- 
rameter of the layer compound is obtained as ol = 
5.6691 (4) A, which is the same as the nominal value of 



5 



TABLE I: Hybrid reciprocal lattice points observed in Figs. 3(b), 3(c), and 3(d), as indicated by HRLP numbers and h'k'l' 
indexes. Qxy values are obtained by using Aa* — Ab* = —5.15 x 10~*A~^ in equations (jU and (|6}. P' = Ms.l + Nl,s is the 
hybrid diffraction vector, and ^ is the cosine director of the secondary beam, equation ((5)l . A few third-order hybrid diffractions 
are also visible: HRLPs 5, 8, and 9. Exact azimuthal ifih'k'e' angles where each HRLP is excited have been calculated as 
described in Appendix A. The crystal azimuthal position during data collection is given at the if column: values with respect 
to the [110] in-plane reference direction. 5 is the diffracted beam takeoff angle as defined in equation Elliptical-shape 
orientation angles of the HRLPs in the reciprocal-space maps. Fig. 3, are given by xe (experimental values), while XG, and 
Xm are theoretical values according to two distinct hypotheses, as explained in the text. Angular values are given in degrees. 
HRLP h'k'e' Q^y (A-i) P' C/ICI Vh'kW V S XE (±2) XG XM 



1 


311 


-9.50 X 10 


-4 


313s + Silt 


+ 1 


27.7043 


27.6 


+0.116 


+75 


+47.6 


+75.2 


2 


311 


+9.50 X 10" 


-4 


31Ii +3l3s 


-1 


27.8258 


27.6 


-0.116 


-75 


-47.6 


-75.2 


3 


200 


-2.92 X 10" 


-4 


202s + 200i. 


+ 1 


28.5022 


27.6 


+0.087 


+90 


0.0 


+90.0 


4 


200 


+2.92 X 10" 


-4 


200l + 202s 


-1 


28.5928 


27.6 


-0.087 


-90 


0.0 


-90.0 


5 


002 






202s + 002i. + 202s 


+ 1 


28.5865 


27.6 










6 


111 


-3.65 X 10 


-4 


113s + llli 


+ 1 


30.1328 


30.0 


-0.056 


+58 


+70.1 


+58.5 


7 


111 


+3.65 X 10" 


-4 


llli, + ll3s 


-1 


30.0750 


30.0 


+0.056 


-58 


-70.1 


-58.5 


8 


002 






llli, +004s + llli. 


-1 


30.0369 


30.0 










9 


004 






Il3s +004i. + ll3s 


+ 1 


29.8584 


30.0 










10 


331 


-15.4 X 10" 


-4 


333s + 33li. 


+ 1 


44.6644 


45.0 


-0.136 


+80 


+33.4 


+77.7 


10 


331 


-15.4 X 10" 


-4 


333s + 331i. 


+ 1 


45.3356 


45.0 


+0.136 


+80 


+33.4 


+77.7 


11 


331 


+ 15.4 X 10" 


-4 


33li + 333s 


-1 


44.5215 


45.0 


+0.136 


-80 


-33.4 


-77.7 


11 


331 


+ 15.4 X 10" 


-4 


33li, + 333s 


-1 


45.4785 


45.0 


-0.136 


-80 


-33.4 


-77.7 



the ZnSe compound, "while the epilayer undergoes an ex- 
pansive in-plane strain e|| = 1.36(54) x 10~^. 

HRLPs in Fig. 3(d), numbers 10 and 11, have the 
largest observed Qxy component, scattering within a the- 
oretical takeoff angle S = ±0.14°, equation ([7]), but still 
accepted by the analyzer system of the diffractometer. 
In fact, the minimum required acceptance of the ana- 
lyzer system should be, in general, equal or larger than 
the axial divergence since the experimental takeoff an- 
gles are increased by an amount corresponding to the 
difference bet'ween iph'k'i' and if values in Table 1. For 
instance, HRLPs 3 and 4 -would not appear on the map 
in Fig. 3(b) if the axial acceptance "was belo"w 1.08°. 

HRLPs "with null Qxy component, such as those num- 
bered 5, 8, and 9 in Figs. 3(b) and 3(c), can not be related 
to second-order sequences of reflections. 00^' hybrid re- 
flections do not occur for partially or fully relaxed epi- 
layers on (001) substrates because Aa* and A&* in equa- 
tion ^ have non-null values. Sequences of third-order, 
such as IIIl + 004s -I- 111^ and 113s + 004^ -I- lT3s 
can explain the 002 and 004 type of HRLP even in re- 
laxed layers. Note that "when carrying out longitudinal 
scans, i.e. cu : 26 scans, HRLPs along the primary re- 
flection axis generate extra features even in samples "with 
a single epilayer. Moreover, in ordinary rocking curves, 
hybrid intensity contributions could be misinterpreted as 
due to sublayers or other sort of structural change, such 
as composition grading. 

T"wo hypotheses have been verified to explain the elon- 
gated shapes of the experimental HRLPs in Fig. 3. The 
strain grading hypothesis should, according to equation 
(IHl) , reduce the absolute value of xg as the in-plane com- 
ponent of P' increases. The expected values of xg, calcu- 
lated using 1/ = 0.5, are given in Table 1. Although elon- 



gations happen in the same sense as the observed ones, 
their behavior as the h'k' indexes increase is the oppo- 
site of that expected. For instance, XG = +70.1°, -1-47.6°, 
and -1-33.4°, are the expected values for the HRLPs num- 
bered 6, 1, and 10, "while the observed values are instead 
XE = +58(2)°, +75(2)°, and +80(2)°. 

A successful explanation for the elliptical shapes of the 
HRLPs comes from mosaicity in the epilayer. Its 002 
reciprocal-lattice point in Fig. 3 (spot labeled L) has a 
fuU width at half-maximum (FWHM) of 17.5 x lO-^A"! 
along the Qxy coordinate, yielding an out-plane mosaic 
width of r]± = 0.024°. Fig. 4(a) shows the computed 
directions of maximum projection of the AH vector onto 
the QzQxy maps, described by the angle xm, equation 
PT|) . as a function of mosaicity-anisotropy ratio r]^\/r]±. 
A very good match to the experimental values is obtained 
for a ratio of 0.67(3), as can be seen by comparing xm 
and in Table 1, and then, 77|| = 0.016(1)°. Projection 
shapes for this ratio value is given in Fig. 4(b). In the- 
ory (Fig. 4a), hybrid reflections of indexes 111 and 111, 
HRLPs 6 and 7, have better sensitivity to the in-plane 
mosaicity, and they are used to estimate the ±0.03 error 
bar of the found ratio value. 



IV. CONCLUSIONS 

Hybrid reciprocal lattices for x-ray diffraction in epi- 
layer/substrate systems exist. Without a systematic de- 
scription, as the one provided in this work, accidental 
excitement of hybrid reflections in most X-ray diffrac- 
tion techniques can jeopardize data analysis. On the 
other hand, understanding their properties leads to three- 
dimensional information of the layer structure even in 
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FIG. 4: (a) Orientation angle xm of the HRLPs 1 and 2 (dashed Une), 6 and 7 (soHd hne), and 10 and 11 (dot-dashed hne), 
Table 1, as a function of the anisotropy ratio ?7||/?7± between the in-plane, 77||, and out-plane, ri±, mosaic widths. Values were 
computed according to equation pi[) . Experimental values (squares) from the reciprocal-space maps in Fig. 3 are also shown, 
(b) Theoretical shapes and orientations for 67% of anisotropy in the epilayer mosaicity, numbered and sized according to Fig. 3. 
Obtained xm values are compared to the experimental ones in Table 1. 



symmetrical high-angle difFraction geometries. In the 
case analyzed here, from reciprocal-space maps of a single 
symmetric reflection we were able to determine parallel 
and perpendicular lattice mismatches, the stress state of 
the epilayer, the absence of strain grading, and the spa- 
tial misorientation of mosaic grains in the layer. 



APPENDIX A: BRAGG CONDITION FOR 
HYBRID REFLECTIONS 

Let us initially consider a given three-beam diffraction 
of the substrate lattice in which P = M + N. In the 
reference system {x, y, z) where i = P/|P|, x is an 
arbitrarily chosen in-plane direction perpendicular to P, 
and y — z X X, the incident beam wavevector 

k — — A^"'^(cos w cos ip X -\- cos (jj sin <^y + sin w z) (Al) 

fulfills the two Bragg conditions in equations ([T]) and ^ 
when Lij = ujQ and ip = Lp^. Hence, wq is the Bragg angle 
of the primary reflection, diffraction vector f, and <po 
can be calculated by the expression 



cos((y9o — a) 



\\M.\I2 — sincjQ COS7 



(A2) 



cos loq sm 7 

given that M — |Af |(sin7COSQ!£-|-sin7sinay-|-cos7i). 

The problem is how to calculate the exact incidence, 
Lo — ~\- Aw, and azimuthal, = 939 + A(^, angles of 
a hybrid reflection whose effective diffraction vector is 
P' = Ms.L + Nl.s- Since P' is not parallel to the z 
axis, equation (jA2p is no longer valid for either Ms or 
Ml diffraction vectors of the secondary reflection. To 
solve this problem we first wrote (Morelhao, 2002) 







dk 
du 



Alo ■ 



dk 

dip 



Aip^ko + k^Au + k^Aip (A3) 



and then equations ([2]) and ([3]) were used to build the 
system of linear equations 



M k^ M 



■N 



k^-N 





Alo 






Atp 





(M/2 + fco) • M 
{N/2 + M + ko)-N 
(A4) 

where M and N stand for Mg ^ and AT^ respectively. 

The iph'k'i' values in Table 1 were calculated by solving 
the above equations. For instance, consider the 111 
HRLP, hne 7 in Table 1. Since the 002 GaAs reflection is 
the primary one and the [110] crystallographic direction 
was chosen as reference for the tp rotation, 

i= [1,1,0]/V2, y = hl,l,0]/V2, and [0,0,1]. 

The corresponding substrate three-beam diffraction oc- 
curs at 

Wo = 15.8132° and ipo = 30.0442°, 



as obtained from equation (jA2p since 7 = 125.2644°, a = 
90°, and |M| = V3/ao. The uq and ipo angles provide 
the wavevector fco from equation (jAip , as well as k^^ and 
fcu,. By replacing in equation (IA4I) 



M 



[-1/a, 1/5,-1/c] and AT ^ Ns ^ [l,-l,3]/ao, 

with a = 5.6699A, c = 5.6676A, and oo = 5.6534A, 
we have 

Auj = 0.0794° and Aip = 0.0308°, 

which leads to iph'k'i' — fo ^ ^'P — 30.0750° in Table 1. 

An alternative approach to calculate the Qz and Qxy 
coordinates of the HRLPs, equation ([S]), is available af- 
ter determining Alo and Aip. Since fc is known, we also 
known k' = P' + k. 



fc.T 



- (cos ip X + sin ip y) and s = (sin ip x ~ cos ip y). 



7 



By defining 

^off = i^' ' ^^'^ Kn = k' k'offi 

we have 

w' = arccos(fc^y ■Knl\Kn\). 

and hence 

= (sinw' + sincj)/A - |P| 
Qi^j/ = [cos uj' — cos uj) I \. 

This approach also provides an alternative expression, 

6 — arctan(Afc' • J), 



to calculate the takeoff angle given in the equation 
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